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Abstract. This article summarizes the basic facts and ideas concerning the formation and 
evolution of cataclysmic variables (CVs). It is shown why the formation of CVs must in- 
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velope evolution. A brief discussion of the principles of the long-term evolution of semi- 
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1. Introduction 

Cataclysmic variables (CVs) are short-period 
semi-detached binary systems in which a white 
dwarf (WD) primary acc retes matter f rom a 
low-mass companion star (Warner 1995). CVs 
are intrinsically variable and that on a wide 
range of time scales (from seconds to > 10 6 yr) 
and with a huge range of amplitudes (of up 
to 10 6 and possibly even more). The rich phe- 
nomenology of CV variability which includes, 
among other things, phenomena like flicker- 
ing, dwarf nova and classical nova outbursts, 
can to a large extent be understood as either 
immediate or long-term consequences of the 
mass transfer process. Interesting as all these 
phenomena are, they are of no particular in- 
terest here. Rather, in the following I shall 
concentrate on evolutionary aspects, i.e. on 
the formation and evolution of CVs. Readers 
who are mainly interested in CVs as variable 
stars should ins te ad turn to the m onographs by 
IWarner I (fl995h or lHellier I d2002l) . 
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2. Very basic facts about CVs and 
stellar evolution 

2.1. Generic properties of CVs 

From the perspective of stellar evolution a CV 
is a semi-detached binary in which a WD pri- 
mary of mass Mi accretes from a low-mass 
secondary star of mass M2 which fills its crit- 
ical Roche lobe. From Roche geometry it fol- 
lows that the secondary's radius can be written 
as R2 - a fi(q). Here a is the orbital separation, 
q - Mi IM2 the mass ratio, and f% the fractional 
Roche radius of the donor star. For typical val- 
ues of q found in CVs, i.e. 1 < q < 10, Eqs.® 
or (01, given below, yield 0.2 < fa < 0.4. 

In principle, the mass of the WD compo- 
nent can be anywhere between the lowest pos- 
sible value resulting from stellar evolution (~ 
0.1 5M ) and the Chandrasekhar mass Men ~ 
1.4M©. Observed masses are mostly in the 
range O.5M < M\ < 1M Q . As to the mass 
distribution there are reasons to believe that in- 
trinsically it is not unlike that of single WDs 
which have a mean mass of < M WD >x O.6M . 
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From the observed mass transfer rates one 
can infer that mass transfer in CVs is sta- 
ble. This, in turn, requires that the mass of 
the donor is typically less than that of the 
WD component, i.e. M 2 < M\, or q > 1, 
and thus that the donor is a low-mass star. 
Observations show that in more than 95% of 
all cases the donor star is on the main sequence 
(MS), though not necessarily close to the zero 
age main sequence (ZAMS). In rare cases the 
donor star is either a giant, or a WD of very 
low-mass (M 2 < O.O5M ). 

For later comparison it is useful to keep in 
mind the resulting typical system parameters 
of a CV with a MS donor: 

- total mass: M = M\ + M 2 ~ M 

- orbital separation: a » few R 2 « R Q 

- orbital period: 80 min < P orb < 10 h 

- orbital angular momentum: 

7 orb = G lll M x M 2 (Mj + M 2 T 1/2 a 112 (1) 
* J = G i/2 M e 3 ' 2 Re 1 ' 2 . 

2.2. Evolution of single and binary stars 

In the following I summarize the basic facts 
which characterize single star and binary evo- 
lution, and which are of relevance in the con- 
text of our considerations. These facts are: 

1. Stars grow considerably as they age (by 
factors up to > 10 2 ). Because this growth is 
not strictly monotonic one can distinguish 
distinct evolutionary phases during which 
a star grows. These phases are: 

- central hydrogen burning, i.e. on the 
MS 

- for intermediate mass and massive stars 
(M > 2.2M ) the post-MS evolution 
towards He-ignition including the evo- 
lution through the Hertzsprung gap 

- for low-mass stars (M < 2.2M ) evolu- 
tion on the first giant branch up to the 
He-flash 

- for low and intermediate mass stars 
(M < 1OM ) evolution on the asymp- 
totic giant branch (AGB) 

2. The more massive a star, the faster it 
ages. Stars on the main sequence obey a 



mass luminosity relation. On the upper MS 
(1M < M < 10 M ) the luminosity L 
scales roughly as L oc M 3 5 . Hence the nu- 
clear time scale is r nuc oc M/L oc M~ 2 5 . 
The immediate consequence of this is that 
of two stars with the same age (as in a bi- 
nary) but different mass the more massive 
star grows faster, i.e. is the bigger of the 
two. 

3. In a binary the presence of a compan- 
ion limits the size up to which a star can 
grow (Roche limit) without losing mass to 
its companion. The maximum radii corre- 
sponding to the Roche limit are the criti- 
cal Roche radii R\$_ - a f(q) and R 2 $_ = 
a f 2 (q) for respectively the primary and the 
secondary ; where f\(q) = f 2 (l/ q), and ac- 
cordin g to lPacz vnski I d 1 97 ll) and lEggleton I 
419831) for 1 < q < 10 

f 2 (q) * 0.462(1 + qy U \q > 1.25 (2) 
f(q) * 0.38 + 0.2 log q*q 0A5 f 2 (q). (3) 

As a consequence of this, stellar evolu- 
tion in a binary of not too large an orbital 
separation a results sooner or later in the 
formation of a so-called semi-detached bi- 
nary in which the more massive component 
reaches its Roche limit first and starts trans- 
ferring mass to its companion. 

2.3. Prerequisites for white dwarf 
formation 

WDs are the end product of the evolution 
of stars of low and intermediate initial mass. 
Thereby the chemical composition of a WD re- 
flects the evolutionary state of the star when 
it loses its hydrogen-rich envelope. Depending 
on when this happens along the evolution the 
result is either a WD consisting mainly of he- 
lium (He-WD), of carbon and oxygen (CO- 
WD), or oxygen and neon (ONe-WD). 

- He-WDs result from the complete loss of 
the hydrogen-rich envelope of a low-mass 
star (with an initial mass Mi < 2.2M ) 
on the first giant branch, i.e. before reach- 
ing the He-flash. Accordingly, the mass 
of He-WDs is in the range 0.15 Af© < 
Mne-WD < Af He -Fb where M He -Fi ~ 0.45 - 
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O.5OM is the mass of the He core at 
the onset of the He-flash. Because wind 
mass loss of single stars on the first gi- 
ant branch is not strong enough for com- 
plete envelope loss isolated He-WDs are 
not formed. However, they can result from 
mass transfer in a close binary (se e e.g. 
iKippenhahn. Kohl & Weigertl dl967l) 1. 

- CO-WDs result from the complete loss 
of the hydrogen-rich envelope of interme- 
diate mass stars on the AGB, i.e. before 
the onset of carbon burning. For single 
stars this happens if the initial mass is 
Mi < 6 - 8M . In binary stars this can 
happen for initial masses up to ~ 1OM . 
Accordingly, the resulting WD masses are 
in the range M He -Fi 5 M C0 -wd < M C -ig n , 
where Mc-i gn ~ 1.1 M Q is the core mass at 
the onset of carbon ignition. 

- ONe-WDs originate from stars which un- 
dergo off-center carbon ignition and sub- 
sequent envelope loss during the so-called 
super-AGB phase. For single stars this is 
possible for initial masses in the range 
9M < M; < 1OM , whereas in bina- 
ries the mass range is 9M < M\ < 
12M (see e.g. Gil-Pons & Garcfa-Berro I 
d200ll) . lGil-Pons et al. I d2003h ). The result- 



ing WDs have masses in the range 1 . 1 M G 



M ON , 



e-WD 



< 1.38M F 



In the context of our considerations, one of 
the most important properties of stars which 
have a degenerate core of mass M c is that they 
obey by and large a core mass-luminosity re- 
lation _£(M C ), and to the extent that these stars 
have a sufficiently massive hydrogen-rich en- 
velope and thus are close to the Hayashi-line, 
also a core ma s s-radiu s relation 7?(Mc) (se e 
e.g. iPaczvnskil dl970l). IKippenhahn I dl98ll) . 
IJoss. R appaport & Lewis I (1198^ ^). This rela- 
tion shows that the radius of such a star is a 
steeply increasing function of core mass and 
that, in particular, AGB stars and stars on the 
super-AGB are very large with radii of up to 
~ 1O 3 R . In other words: the formation of a 
WD requires a lot of space, the more massive 
the WD the more space. This is not a problem 
for single stars. But in a binary, as a conse- 
quence of the Roche limit, the orbital separa- 



tion a sets an upper limit to the mass of the 
WD that can be formed: M WD < R. 1 (af(q)). 

2.4. Single star evolution versus binary 
star evolution 

The task of calculating the structure and evolu- 
tion of a single star consists of solving a well- 
known set of differential equations with appro- 
priat e boundary conditions and ini tial values 
(e.g. lKippenhahn & Weigertl d 19901) ). 

For calculating the evolution of a binary 
system (or of one of its components) the task 
is in principle the same as for single stars. The 
difference is that in a binary one has an addi- 
tional boundary condition which derives from 
the presence of the companion star, i.e. from 
the Roche limit. 

Consider for simplicity a system consist- 
ing of a "real " star, say the primary, and a 
point mass secondary. The simplest boundary 
condition that one could impose in this case is 
that Ri < afi(q). A more realistic approach 
would take into account that the surface of a 
star is not arbitrarily sharp, but rather is char- 
acterized by a finite scale height H <K R over 
which pressure, density etc. drop off, by ex- 
pressing the mass loss rate -M\ as an explicit 
function of b inary and stellar parameters (e.g. 
iRitterl dl988l) ). What one finds is that -M x is 
a steeply increasing function of {R\ - Rir)/H 
and that the primary suffers significant mass 
loss as R\ — > Si r. 

The real problem when dealing with mass 
transfer consists of answering two questions: 
1.) Where does the mass lost from the donor 
go? and 2.) How much angular momentum 
does it take with it? On the formal level this 
can be dealt with as follows: let us assume that 
a fraction tj of the transferred mass is accreted 
by the secondary, i.e. 



M 2 - -Tj Mi . 



(4) 



Accordingly, the mass loss rate from the sys- 
tem is M, = (1 — rj) Mi. The angular momen- 
tum loss rate associated with this mass loss can 
be written as 



j olh = vMJ orb /M, 



(5) 
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where v is a dimensionless factor measuring 
the angular momentum leaving the system. 

What is known about the two parameters rj 
and v? In general not more than 1.) < t] < 1, 
and 2.) v > 0. Otherwise rj and v are (al- 
most) free functions of the problem. Therefore, 
calculating binary evolution involves (at least) 
two almost free functions. And the art of per- 
forming such calculations very often consists 
of making creative use of this freedom! 

2.5. Generic properties ofCV 
progenitors 

We are now in a position to define the neces- 
sary criteria which a binary consisting initially 
of two ZAMS stars of mass My and M24 has 
to meet in order to later become a C V which, at 
the onset of mass transfer, i.e. at the beginning 
of its life as a CV, consists of a WD of mass 
Mwd and a donor star of mass M2. 

1 . M\i has to have sufficient mass to allow for 
the formation of a WD of mass Mwd- 

In theoretical calculations of the evolution 
of single stars with a fixed set of phys- 
ical assumptions (such as initial chemi- 
cal composition, equation of state, opaci- 
ties, nuclear reaction rates, convection the- 
ory, wind mass loss, etc.) there is a one to 
one relation between the initial mass Mi 
and the mass Mf of the white dwarf pro- 
duced. This relation is known as the ini- 
tial mass-final mass relation, i.e. Mwd — 
Mf(M;). And, within the observational un- 
certainties, there is also ample observa- 
tion al evidence for t his M\ -Mf -relation (see 
e.g. ISalaris et al. I ((2008) and references 
therein). 

In binary evolution things are a little differ- 
ent: because mass transfer sets a premature 
end to the donor's nuclear evolution the 
mass of the resulting white dwarf is smaller 
than what single star evolution of the pri- 
mary would yield, i.e. Mwd < MffMy). 
In other words: for the formation of a WD 
of mass Mwd the necessary condition is 
M u > Mr 1 (M WD ). 

2. Because of the core mass-radius relation 
!R(M C ) which holds for the giant primary 



when it reaches its Roche limit, the ini- 
tial separation of the binary must be a\ = 
•R(MwD)//i(<7iX where q\ is the initial mass 
ratio. For this estimate of a\ we have im- 
plicitely assumed that after the onset of (the 
first) mass transfer Mwd = const. 
3. Finally for the secondary's mass we as- 
sume M2,i = M2. A justification for this 
will be given below. 

Now, let us take typical parameters for a CV, 
say M W d ~ 1M Q and M 2 < 1M , in order to 
see where this leads us: with Mwd ~ 1M it 
follows from the Mj-M f -relation that My > 
5M Q , hence Mi > 6M Q , and from the core 
mass-radius relation 'R(Mwd) ~ 1O 3 7? , and 
with /i(<7i) ~ 0.5, a; ~ 2 1O 3 7? . Therefore, the 
initial orbital angular momentum of the binary 
is 

« io 2 / . 

Comparing now the total mass and orbital an- 
gular momentum of a CV (cf. Sect. 12. U with 
the corresponding values of its progenitor sys- 
tem we find that Mi/Mcv ~ 5-10 and 
J ib,i/Jcv ~ 10 2 . In other words: the formation 
of a CV invokes a binary evolution in which the 
progenitor system has to lose ~ 80% - 90% of 
its initial mass and up to ~ 99% of its initial 
orbital angular momentum (lRitterlll976l) . and 
that after the onset of mass transfer from the 
primary. 



3. Mass transfer and its 
consequences 

Since the primary of a CV progenitor does 
not stop growing when approaching its Roche 
limit, onset of mass transfer is unavoidable. 
And, because the subsequent formation of a 
CV involves huge losses of mass and orbital 
angular momentum from the binary system, it 
is necessary to examine the consequences of 
mass transfer for the ensuing evolution in more 
detail. 
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3.1. Stability of mass transfer 

A detailed discussion of the stability of mass 
transfer is rather complex and beyond the 
scope of this article. For th is t he read e r is re - 
ferred to e.g. iRitter I dl988l) or lRitterl dl996l) . 
Here I shall keep the presentation as simple as 
possible. 

Let us assume for the moment that the pri- 
mary star has a sharp outer boundary and that 
it has just reached its Roche limit, i.e. that 
R\ = Ri,r- What happens if at that moment, 
which we denote by to, a small amount of mass 
6m is taken away from the primary and trans- 
ferred to the secondary, i.e. if Mi — » M\ — 5m 
and Mi — > M2 + 5ml As a consequence of 
this small mass transfer, not only the mass ra- 
tio q and the critical Roche radii R\ n and /?2,r 
will change but also the stellar radii R\ and R2- 
Let us for the moment treat the secondary as a 
point mass. Then we have to deal only with the 
radii R[(t > to) and Ri^(t > to). Thereby, three 
different situations can arise: 

1. Ri(t > to) < ^i,r(? > to): In this case mass 
transfer is stable, because after a small 
mass transfer 5m the donor underfills its 
critical Roche volume and mass transfer 
stops. 

2. R\(t > to) > R\,R(t > to): In this case mass 
transfer is unstable, because if R\(t > to) — 
Ri,r(1 > h)) > even more mass flows over. 

3. Ri(t > to) - /?i,r(? > fo): In this case mass 
transfer is marginally stable. 

In order to decide which of the three above 
cases arises we must know how Ri and /?ir re- 
act to mass transfer. For all practical purposes 
Z?i_r adjusts instantaneously (actually on the 
orbital time scale) to changes in Mi, M2 and 
/ or b. Although, in principle, calculating Ri r is 
straighforward, for this it is still necessary to 
precisely specify where the transferred mass 
goes and, if the system loses mass, how much 
angular momentum it takes with it, i.e. one has 
to specify the parameters 77 and v. The change 
of Si r is conveniently expressed in terms of 
the mass radius exponent 



&.1 



(dlnR 



l.R 



<91nM, 



(7) 



where the subscript * is a reminder that for its 
calculation 77 and v need to be specified. 

On the other hand, the reaction of the 
donor's radius Ri to mass loss is more compli- 
cated: besides hydrostatic equilibrium which 
readjusts on the orbital time scale, mass loss 
disturbes also the thermal equilibrium of a star. 
Therefore, its reaction depends on the ratio of 
the mass loss time scale tm to the time scale 
T t h on which the star can readjust to thermal 
equilibrium. If tm/t^ <k 1 the star reacts es- 
sentially adiabatically, and the radius change is 
expressed in terms of the adiabatic mass radius 
exponent 



£ad,l - 



dlnR, 



din Mi 



(8) 



l/ad 



If, on the other hand, mass loss is very slow, i.e. 
tm/t± » 1, the star has time to adjust to near 
thermal equilibrium in which case the radius 
change is expressed by the thermal equilibrium 
mass radius exponent 



dlnR, 



dlnM, 



(9) 



i/th 



Accordingly, there are two criteria for the sta- 
blity of mass transfer: 

1 . Mass transfer is adiabatically stable if 

<Tad,l - <Tr,1 > (10) 

2. Mass transfer is thermally stable if 



£h,i - &,i > . 



(11) 



What does all that mean for the CV progeni- 
tor system at the onset of mass transfer? In or- 
der to tell one has to know the values of £r_i, 
£ a d,i, and fth,i- Because M14 > M24 one in- 
variably finds that (rj > even in the most 
favourable case where no orbital angular mo- 
mentum is lost. The values of £,d, 1 , and ^th, 1 , on 
the other hand, depend on the internal structure 
of the star in question. In our case the donor 
is a star with a degenerate core and a deep 
outer convective envelope. For such stars one 
typically finds -1/3 < £,rt ^ and £ t h ^ 
(Hiell ming & Webbink|[l987l) . Taken together 
this means that mass transfer in such a system 
is adiabatically and thermally unstable. And 
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as a consequence of the adiabatic instability 
mass transfer quickly accelerates to the point 
where the mass transfer rate reaches values of 
order of -Mi jac | ~ Mi/r conv ~ M yr _1 , where 
Tconv ~ yr is the convective turno ver time scale 
dPaczvhski & Sienkiewicziri972l) . 

3.2. Fast accretion onto a main 
sequence star 

So far we have treated the MS secondary 
as a point mass. Whereas before the on- 
set of mass transfer this is an adequate 
approximation, this is not always true 
afterwards. Numerical calculat ions (e.g . 
Kippenhahn & Mever-Hofmeister I d 19771) . 
Neo et al. I (11977b ') show that the low-mass 
secondary, exposed to the prodigious mass 
inflow rates associated with the adiabatic 
mass transfer instability, starts expanding 
rapidly to giant dimensions. The reason for 
this behaviour is that the thermal time scale of 
the accreted envelope around the secondary is 
much longer than the mass accumulation time. 
As a consequence, the accreted matter can not 
cool efficiently and, therefore, forms a deep 
and very extended convective envelope of high 
entropy material around the secondary. The 
star thus attains a structure similar to that of 
a giant/AGB star which, however, derives its 
luminosity mainly from accretion rather than 
from nuclear burning. 

3.3. Formation of a common envelope 

The situation of a CV progenitor at the on- 
set of mass transfer can now be character- 
ized as follows: because mass transfer occurs 
from the more massive star, the orbital sep- 
aration a as well as the critical Roche radii 
Z?i_r and /?2,r shrink. At the same time, the 
mass losing donor star has the tendency to ex- 
pand (negative £ a d and But forced by dy- 
namical constraints to essentially follow /?i r 
the donor must lose mass at rates approach- 
ing ~ M yr~'. And the secondary, in turn, ex- 
posed to such enormous accretion rates, reacts 
by rapid expansion. The consequence of all this 
is that within a very short time after the onset 



of mass transfer the system evolves into deep 
contact. An attempt to model this very com- 
plicate d process has been made by IWebbinkl 
(Il979l) . Accordingly, the immediate result of 
this evolution can then be roughly chracterized 
as follows: A binary system consisting of the 
primary's core (the future WD) of mass M c and 
the original secondary of mass M2;, finds itself 
deeply immersed in a common envelope (CE) 
of mass Mce = - M c and a size which 
must be of order of or even larger than the ra- 
dius given by the core mass-radius relation, i.e. 
R CE > 7?(M C ). 

4. Common envelope evolution and 
CV formation 

Common envelope evolution is the name of a 
generic process which arises as a consequence 
of dynamical time scale mass transfer and as 
a result of which a detached short-period bi- 
nary is formed in which one of its compo- 
nents is the core of the former primary (in 
our case a pre- WD). Because of its importance 
for the formation of all sorts of compact bi- 
naries the subject has generated a vast liter- 
ature. For lack of space I am unable to give 
a detailed review here. Rather I shall concen- 
trate on sketching a few key aspects of this pro- 
cess and for mo re details refer the reader t o re- 
cent reviews bylTaam & Sandquist I d2000h and 
IWebbinkld2008h . 

4.1. The Darwin instability 

Let us now consider the following idealized 
situation: a binary consisting of the original 
primary's core of mass M c and the secondary 
of mass M2 with orbital separation a and or- 
bital frequency a>„ is embeded in an envelope 
of mass Me, radius Re, moment of inertia Ie 
which is in solid body rotation with an angu- 
lar frequency Qe- If <w« > Qe tidal interac- 
tion and friction between the binary and en- 
velope lead to energy dissipation and angular 
momentum transport from the binary to the 
envelope with j„ = -Je < 0. As a conse- 
quence, the envelope, initially rotating slower 
than the binary, is spun up. But according to 
Kepler's third law also the binary's orbital fre- 
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quency increases due to the loss of orbital an- 
gular momentum. The question of interest is 
thus whether through this spin-up the differ- 
ence fc>„* - Q E increases or decreases. 

If 6J«» - Q E > and co„ - Q E < the 
envelope is synchronized, i.e. Q E — > 

If, on the other hand, a>„ - Q E > and 
d>„ - Q E > 0, runaway friction results, and the 
binary spirals in. The condition for this to hap- 
pen is easily derived: the binary's orbital angu- 
lar momentum is 

M C M 2 _ 1/3 



G 2/3 



(m c + M 2 y /3 



= h»OJ» 



(13) 



(14) 



where 

M C M 2 2 

h* a 

M c + M 2 

is the orbital moment of inertia. The envelope's 
spin angular momentum is 

/e = / B Ob. (15) 

With (TT~2T > and dT3T > angular momentum conser- 
vation, i.e. j„ + 7 E = 0, yields 



Q F 



w« 1 



3 h 



(16) 



From (TT&b it is seen that the envelope can be 
synchronized only if / E < l/3/»». If, on the 
other hand, 



4.2. Common envelope evolution 

Despite decades of heroic efforts to model 
comm on envelope evolution, f or a review see 
e.g. lTaam& Sandquistl d2000l) , to this day it 
has not yet been possible to follow such an evo- 
lution from its beginning to its end with really 
adequate numerical computations. Therefore, 
it is still not possible for a given set of initial 
parameters to reliably predict the outcome of 
common envelope evolution. The expectation 
is that in many, but not necessarily all cases 
the frictional energy release will unbind the 
CE and leave a close binary consisting of the 
former primary's degenerate core and the sec- 
ondary. 

Clearly the ejection of the CE requires the 
release of the envelope's binding energy in a 
sufficiently short time, i.e. that the time scale of 
the spiral-in is short. However, there are limits 
to how short the spiral-in can be. From sim- 
plified on e-dimensional hydrostatic mo d el cal- 
culations iMever & Mever-Hofmeisterl d 19791) 
found that there is a negative feedback between 
the frictional energy release and the resulting 
radiation pressure. An estimate of the duration 
of the spiral-in is obtained from the argument 
that because of this feedback the frictional lu- 
minosity Lfrict can not exceed the Eddington lu- 
minosity 



/e>3^ (17) L Edd 

the envelope cannot be synchronized and 
spiral-in of the binary is unavoidable. The im- 
possibility of synchronizing the envelope re- 
sults from a variant of an instability which is 
actuall y long known: discovered by iDarwin I 
(Il879h . though in a different context, it is com- 
monly called Darwin instability. 

Whether the Darwin instability is of rele- 
vance for our problem, i.e. whether the crite- 
rion ( fT7T > is met with the formation of a CE 
after the onset of adiabatically unstable mass 
transfer, needs of course first to be checked. 
Since adequate model calculations of the for- 
mation of a CE are still not feasible, simple 
estimates must do. And these indicate indeed 
that for typical parameters of CV progenitor 
systems the forming CE systems are Darwin 
unstable. 



AnGcM 



(18) 



by much. Here K es is the electron scatter- 
ing opacity. The evolution of the binary with 
masses M c and M 2 from an initial separation 
a\ to a final separation af «: a\ releases the or- 
bital binding energy 



GM C M 2 
2cif 



(19) 



This yields a rough estimate of the spiraling-in 
time scale 



1"CE 



ifrict 

> 400yr 



^Edd 

M C M 2 _ 

(M C +M 2 )M Q R e 



(20) 
(21) 



Thus for the typical parameters of a CV (see 
Sect l2.U tqe is very short, so short indeed that 
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the secondary star has no time to accrete a sig- 
nificant amount of mass during the CE phase 
dHiellming & T aam 1 1199 lh . This is the a pos- 
teriori justification for our assumption in Sect. 
1231 that M 2 ,i = M 2 , f . 

Because of the short duration of CE evolu- 
tion chances of observing a binary system dur- 
ing this phase are extremely small, apart from 
the fact that it is not even quite clear what to 
look for. Worse, the spiraling-in binary is hid- 
den from view as long as it is inside the CE. In 
view of our limited theoretical understanding 
of CE evolution in general and the ejection of 
the CE in particular, and the fact that this pro- 
cess is virtually unobservable, one has to ask 
why we can be sure that CE evolution really 
happens as described above. Beyond all the un- 
certainties, the concept of CE evolution does 
make at least one prediction that is testable: 
at the end of the CE process, if the envelope 
is ejected, we expect a binary inside the now 
more or less transparent envelope. And in this 
binary the primary's degenerate core emerges 
as a very hot pre-WD which, in turn, ionizes 
the surrounding gas, thereby transforming the 
ejected CE into a planetary nebula. The con- 
cept of CE evolution thus implies the exis- 
tence of planetary nebulae with short-period 
binary central stars. And indeed, such objects 
are obeserved: currently we know of ~ 20 
short-period bin ary central stars of planetary 
nebula e ("s ee e.g. | De Marco. Hillwig & S mith 
(120081) and lRitter & Kolb I (l2003h V 

4.3. Formal treatment of the CE phase 

CE evolution, if it ends with the ejection of 
the CE, transforms a binary with initial param- 
eters (Mi i, M2 ; i, a;) to one with final param- 
eters (Mi_f, M 2) f, ai). With current theory it is 
not possible to precisely link these two sets of 
parameters. Therefore, in evolutionary studies 
and population synt hesis cal c ulation s of com- 
pact b in aries (e.g. Ide Koofl dl990h . Ide Kool I 
1992), de Kool & Ritter (1993). iPolitanol 



1996), Pol itanol (l200l . lPohtano I d 2007t) ). CE 



evolution is usual ly dealt with by means of 
a simple estimate (Webbink 1984) which de- 
rives from the assumption that a fraction »ce ^ 
1 of the binary's binding energy which is re- 



leased in the spiraling-in process, AE^,**, is 
used to unbind the CE. 

Using M u = M c ,i = M c , M 2 ,t = M 2ji = M 2 
we have 



GM C M 2 / 1 



1 

Of 



(22) 



On the other hand, the binding energy of the 
CE can be written as 



iB.CE 



G Mi,j M CE 



(23) 



where Mqe = M\j - M c is the mass and 
Ri,i = fli/i(<?i) the radius of the CE, and A 
a dimensionless factor which can be deter- 
mined from stellar structure calculations pro- 
vided one knows exactly where the mass cut 
between core and envelope is. Unfortunately it 
turn s out that A depends rather sensitively on 
this dTauris & Dewi 1 1200 ll) . The CE criterion, 
namely that 



£b,ce = a A£ B ,*» 
is then equvalent to 



2 Mi j Mce 



M,i 



a C E^M c M 2 fi(qi) M c 



(24) 



(25) 



Eq. d25l l provides the formal link between the 
pre-CE and the post-CE binary parameters. As 
can be seen from Eq. d25l ) when dealing with 
CE evolution in this way one introduces es- 
sentially one free parameter, namely acE A (per 
CE phase). Since we do not have any a priori 
knowledge about q-ce and since also A is not re- 
ally well known, the degree of uncertainty in- 
troduced via acE A is quite considerable. 

Several recent investgations of binary evo- 
lution involving CE evolution have come to 
the conclusion that the energy criterion d24l 
is not always adequate and that in addition to 
the orbital binding energy possibly also other 
sources of energy such as the ionization energy 
have to be taken into account. For a compre- 
hensiv e discussion of this point see Webbink 
(200§. 
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4.4. Evolution of post-common envelope 
binaries 

The ejection of the CE leaves a detached short- 
period binary inside a planetary nebula which 
is excited by the hot pre- WD component. Once 
the planetary nebula disappears, either because 
it dissolves or because of lack of ionizing ra- 
diation from the pre-WD, what remains is a 
binary consisting of a WD and an essentially 
unevolved companion. And because the life- 
time of a typical planetary nebula of ~ 10 4 yr is 
much shorter than the lifetime of a typical post- 
CE binary in the detached phase, the intrinsic 
number of detached post-CE systems lacking a 
visible planetary nebula must be vastly larger 
than that of post-CE systems with a planetary 
nebula. And although such systems are intrin- 
sically rather faint (both the WD and its low- 
mass companion are faint), because of their 
rather high space density quite a number of 
such systems are kn own (currently > 50, see 
iRitter & Kolbl (120031) for a compilation). They 
are collectively refrerred to as pre cataclysmic 
binaries, hereafter pre-CVs. 

In the following, we need to discuss two 
questions: 1) how does a detached pre-CV be- 
come semi-detached, i.e. a CV, and 2) whether 
with the onset of mass transfer all pre-CVs re- 
ally become CVs or perhaps follow a totally 
different evolutionary path. 

Since in a detached system the future donor 
star underfills its Roche lobe, mass transfer 
can only be initiated if either the donor star 
grows (as a consequence of nuclear evolution) 
or if the orbital separation shrinks as a con- 
sequence of orbital angular momentum loss 
(AML). Which of the two possibilities is rele- 
vant for a particular binary system depends on 
the ratio of the nuclear time scale 

Tnuc,2 = ( I (26) 



ldln* 2 / nuc 

on which the star grows to the AML time scale 



at 



= -2 



dt 
<91n< 



(27) 



\d\nJ 0lh/ 

on which the orbital separation a shrinks. 

If tj < 2 T nuCi 2 mass transfer is initiated by 
AML, otherwise by nuclear evolution. The typ- 
ical future donor star of a pre-CV is a low-mass 



MS star. Thus T nuCj 2 > 10 9 yr. AML in such bi- 
naries results ei ther from the emission of grav- 
itational waves (Kraft, Math ews & G reenstein 
1962) or from magnetic braking, i.e. a mag- 
netically coupled stellar wind from the tidally 
locked companion. In typical pre-CV sys- 
tems AML is dominated by magnetic brak- 
ing. Unfortunately, for that case there is as yet 
no theory which would allow computation of 
7 01 b from first principles. Again, simple semi- 
empiri cal estimates (e.g. IVerbunt & Zwaanl 
(1981) ) or simplified t heoret ical approaches 
(e.g. iMestel & Spruit! dl987l) ) must do. For 
the typical pre-CV with a low-mass MS com- 
panion, these estimates yield tj ~ 10 8 yr. 
Thus, for such systems mass t ransfer is typi - 
cally initiated via AML (see e. g. [Ritter1(ll986h . 
ISchreiber & Gansicke I d2003l) ). But the sim- 
ple fact that we do observe a number of long- 
period CVs with a giant donor shows that mass 
transfer can also be initiated by nuclear evolu- 
tion of the future donor star. However, the frac- 
tion of pre-CV systems ending up with a gi- 
ant donor is small and, unfortuna tely, strongly 
model-dependent (Ide Kool 111 9921) . 

When the secondary reaches its Roche 
limit and mass transfer sets in stability of mass 
transfer becomes again an issue. Whether mass 
transfer is stable depends on whether the cri- 
teria which we had derived in Sect. 13.11 but 
now applied to the secondary star, are fulfilled. 
Why is this important? Observations and theo- 
retical arguments show that in the vast majority 
of CVs mass transfer is thermally and adiabat- 
ically stable. In other words: only those pre- 
CVs for which £ui,2-ft,2 > and £th,2-&,2 > 
can directly become CVs. What happens to the 
rest? That depends mainly on the evolutionary 
status of the donor and the binary's mass ratio. 
If we distinguish for simplicity MS stars and 
giants as possible donor stars, then the follow- 
ing cases can arise: 

1. MS donor, mass transfer thermally and 
adiabatically stable — > short-period CV 
(^orb ^ 0.5 d) with an unevolved donor. 

2. MS donor, mass transfer adiabatically 
stable but thermally unstable — » ther- 
mal time scale mass transfer, WD with 
stationary hydrogen burning, system 
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appears as a supersoft X-ray sourc e 
(see e.g. van den Heuvel et al. I d 19921) . 
ISchenker et al. I d2002l) )~ — > CV with an 
artificially evolved MS donor. 

3. MS donor, mass transfer adiabatically un- 
stable — » very high mass transfer rates, sec- 
ond common envelope?, coalescence? 

4. giant donor, mass transfer thermally and 
adiabatically stable — > long-period CV 

(^orb > 1 d). 

5. giant donor, mass transfer either thermally 
or adiabatically unstable — > very high mass 
transfer rates, second common envelope?, 
formation of an ultrashort-period detached 
WD+WD binary? 

5. CV evolution 

CV evolution is a complex subject. Yet, be- 
cause of space constraints, here I can only 
present a brief outline of this topic. For read- 
ers wishingtolearn more about it the reviews 
by lKingl dl988t) and lRitterl i 19961) are a good 
starting point. 

5.1. Mass transfer in semi-detached 
binaries 

If mass transfer in a binary is thermally and 
adiabatically stable, as in the majority of CVs, 
no mass transfer occurs unless some external 
force drives it. And in CVs the driving agents 
are the same as in pre-CVs (cf. Sect, \4.4i , 
i.e. AML and nuclear evolution of the donor. 
Furthermore, if mass transfer is stable and the 
strength of the driving changes only on long 
time scales, mass transfer will be essentially 
stationary. In that case the donor's radius R2 
and its Roche radius /?2,r are equal to within 
very few at mospheric scale heights H <K R2 
(Ritterl ll988h . Thus, to a very good accuracy 
we must have R2 = R2,r, or, using R2 = R2,r, 



dlnR2 d In /?2,r 



(28) 



dt dt 

Now, the donor's radius can change because of 
mass loss, nuclear evolution, and thermal read- 
justment. As mentioned earlier (Sect. l3.lT l mass 
loss (if nothing else) drives a star out of ther- 
mal equilibrium. If mass loss were stopped the 



star evolved back towards thermal equilibrium, 
thereby changing it radius initially at a relative 
rate 



d\nR 2 \ 
~dTl 



1 

Tth,2 



(29) 



where r t h,2 is the thermal time scale. Thus the 
rate of change of R2 can be decomposed as fol- 
lows: 



d\nR 2 M 2 1 1 

Z — _ TT bad,2 H 1 

dt M 2 T th ,2 T nu , 2 



(30) 



On the other hand, the donor's Roche radius 
can change because of mass transfer and AML. 
With (|23 we have 



dlnR 2 , R _Mif _ 2_ 
dt ~ M 2 ' 2 Tj 



(31) 



Eqs. (|28l), d30j>, and (ED finally yield the mass 
transfer rate 



-M 2 = 



1 



■i.-!-4U 



£ad,2 - ^R,2 Vth,2 T nuC; 2 Tj ) 



If mass transfer is sufficienly slow such that the 
donor remains close to thermal equilibrium, its 
radius changes according to 



dlnR 2 M 2 1 

= -TT- £h,2 + 



dt M 2 " ' T nuc , 2 
and together with d28l ) and OTl i we can write 
1 / 1 



(33) 



-M 2 = '- I + — 

£th,2 _ £r,2 \T n uc,2 r J 



(34) 



From what has been said so far, it is easily seen 
that for the sign of the mass transfer rate to be 
correct, i.e. for -M2 > 0, the denominator in 
d32l and d34l l must be positive, i.e. that 



£ad,2 - <Tr,2 > 

and 

6h,2 - &,2 > . 



(35) 



(36) 



With Eqs. (f35T > and d36l l we have thus recovered 
the stability criteria for mass transfer. What this 
implies is that for mass transfer to be stationary 
it has also to be stable. 
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5.2. Computing the evolution of a 
semi-detached binary 



Can we use Eqs. (l32l or d34l i for calculat- 
ing the evolution of a semi-detached binary? 
Unfortunately, this is in general not the case. 
The virtue of Eqs. d32l or ( f34T > and the reason 
why we have derived them here is that they 
show clearly how the long-term evolution of 
a semi-detached binary works: mass transfer 
must be stable and be driven by some mech- 
anism. The obvious ones are the growth of the 
donor star due to nuclear evolution or AML 
which shrinks the binary. A less obvious driv- 
ing agent is the growth of the donor star as a 
consequence of thermal relaxation (cf. d32t). 
However, thermal relaxation, itself mainly be- 
ing a consequence of mass loss, cannot main- 
tain mass transfer for times long compared to 
Tth without external driving by one of the other 
mechanisms. 

The reason why we cannot use Eqs. d32l or 
( f34T > for evolutionary computations is that most 
of the quantities appearing in these equations 
are not explicitly known. In partcular, £ a d, 
T nuc , and r t h require knowledge of the complete 
internal structure of the donor star, i.e. nothing 
less than the whole past history of the binary 
system. Worse, even if all that were known, the 
above quantities can only be determined nu- 
merically. Furthermore, computing £r requires 
specification of v and 77 (Eqs. (01 and (0]l). 
Finally, apart from gravitational radiation, the 
AML rate is not well known and in some cases 
only give n as an implicit function of binary pa- 
rameters dMestel& Spruit Ill987h . Even more 
exotic effects such as irradiation of the donor 
star or the accretion disc can strongly affect 
the quantities ap p earing in ( f32T > or (f34-b (see 
e.g.lRitter I (fl9^ . lBuning&R itter d2004 . or 
iRitter] ((2008J) for more). 

Application of Eqs. (f32t or ( f34-b for 
evolutionary computations is therefore 
limited to cases where the donor star can 
either be approximated by a particularly 
simple stellar model, e.g. by a polytrope 
dRappaport. Joss & Webbink II 19821). a bipolv- 



trope (e.g. iRappaport. Verbunt & Joss I 



(119831) . iKolb & Ritter 
stellar structure data 



d!992l) ), or where 
determined before- 



hand from single star evolution can b e use d 
(Webbink. Rap paport & Savonije I dl983l) . 
IRitter I d!999l) V 

In general, such simplifications are unsatis- 
factory. For a more realistic simulation the full 
stellar structure problem must be solved as de- 
scribed in Sect. 12.41 Stellar evolution is an ini- 
tial value problem. So for setting up a simu- 
lation of a CV evolution one has first to de- 
cide at which moment of the evolution to start 
the calculation, e.g. at the onset of mass trans- 
fer from the secondary, and then to specify at 
least the masses of the components and the in- 
ternal structure, i.e. the evolutionary status of 
the donor star, but as the case may be also the 
structure of the accreting WD. Furthermore, 
one has to adopt values or prescriptions for v 
and 77, and finally to decide what to do about 
AML, in particular about magnetic braking, 
i.e. which of the vari ous prescriptions availabl e 
in the lierature (e . g.lVerbunt & Zwaan I dl98ll) . 
iMestel & Spruit I d 19871) 1 to use. When every- 
thing is set up calculating the evolution is in the 
simplest case just a single star evolution for the 
donor star with variable mass where the mass 
loss rate is an eigenvalue of the problem and 
is determined by the additional outer boundary 
condition , e.g. by R2 < R2,r- 

5.3. A sketch of CV evolution 

The orbital period P or b is the only physical 
quantity which is known with some precision 
for a large nu mber of CVs, current ly for over 
700 objects dRitter & Kolb I 120031) . Reliable 
masses, on the other hand, are known, if at 
all, only for a very small minority of CVs. 
Therefore, much of the work on CV evolution 
in the past 30 years has concentrated on un- 
derstanding the observed period distribution of 
CVs. Broadly speaking, this distribution is bi- 
modal with ~ 45% of the objects having peri- 
ods in the range 3 h < P or b < 16 h , another ~ 
45% with 80min < P orb < 3 h , and the remain- 
ing ~ 10% with 2 h < P orb < 3 h . The dearth of 
objects in the period interval 2 h < P or b < 3 h is 
known in the literature as the period gap. 

The maximum period of ~ 16 h is easily un- 
derstood as a consequence of the facts that 1) 
the donor is a MS star, 2) the mass of the WD 
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is Mwd < Mai ~ 1.4M and 3) mass transfer 
must be stable. 

The minimum period of ~ 80 min, in turn, 
is at least qualitatively understood as a conse- 
quence of mass transfer from a hydrogen-rich 
donor which is mainly driven by gravita- 
tional r adiation dP aczvriski & Sienkiewicz 
(Il98ll). |Paczvnski& Sienkiewiczl d 19831) . 
Rapp aport. Joss & WebbinkT (1198^ . 
Because of mass loss, of the order of a 
few 10 _11 M o yr _1 , the donor star becomes 
more and more degenerate when M2 < 0.1 M 
and its structure changes from that of a 
low-mass MS star to that of a brown dwarf. 
Thereby its effective mass radius exponent 
<Teff,2 = din R2/ din M2 changes from ~ 0.8 
on the MS to -1/3. P 01 -b is minimal when 
<Teff,2 = +1/3. Whether mass transfer near 
the period minimum is really driven by 
gravitational radiation only is currently under 
dispute because of the mismatch between 
the corresponding theoretical prediction 
for the minimum period of ~ 70 min and 
the observed v al ue of ~ 7 7 min (see e.g. 
iRenvoize etaTl d2002l) or iBarker & Kolbl 
(2003) for a discussion). 

The period gap is more difficult to ac- 
count for. Over the years a number of differ- 
ent hypotheses have been put forward to ex- 
plain it. For lack of space I cannot review 
them all here. Rathe r I shall concen t rate o n 
the one hypothesis (ISpruit & Ritterl dl983l) . 
iRappaport. Verbunt & Joss I ([1983)) which, in 
my view, still provides the most plausible ex- 
planation for what we see, and which is known 
in the literature as the disrupted (magnetic) 
braking hypothesis. It postulates that as long 
as the donor star has a radiative core "magnetic 
braking" is effective and CV evolution is driven 
by a high AML rate due to "magnetic braking" 
and gravitational radiation, but that, as soon as 
the donor star becomes fully convective, "mag- 
netic braking" becomes ineffective and thus the 
evolution is driven by AML from gravitational 
radiation only. In the following I shall try to 
explain step by step how the gap arises in the 
framework of this hypothesis. 

First it is important to note that the evolu- 
tion of CVs with a MS donor driven by AML 
leads from longer to shorter orbital periods. 



And MS donor stars with a mass < 1M have a 
convective envelope and a radiative core. With 
decreasing mass, i.e. P or b, the mass of the ra- 
diative core shrinks until at a particular mass 
A^2,eonv, i-e. orbital period P or b = P u , the donor 
becomes fully convective. According to the 
above hypothesis, at this point the AML rate 
drops from a high value which is mainly due 
to "magnetic braking" to a small value due to 
graviational radiation only. 

If "magnetic braking" is sufficiently strong, 
then for periods > P u mass loss from the donor 
occurs on a timescale much shorter than its 
thermal time scale. As a result the donor is 
significantly driven out of thermal equilibrium 
and, therefore, oversized compared to its ther- 
mal equilibrium radius, i.e. P2(Porb > Pu) > 
7?2,e, and the faster the mass loss, the larger 
the difference R2 - /?2,e- Suppose now that the 
driving AML rate drops by a large factor on a 
short time scale. What will happen? The donor 
will detach from its Roche lobe because ini- 
tially it will continue losing mass and shrink 
at the same rate as before while its Roche ra- 
dius, because of the reduced AML rate will 
shrink much more slowly. So mass transfer 
stops and the star, being oversized because of 
previous high mass loss but now without mass 
loss contracts towards its thermal equilibrium 
radius /?2,e, and that on its thermal time scale 
which is initially shorter than the time scale 
on which its Roche radius shrinks. Mass trans- 
fer can only resume when the shrinking Roche 
radius reaches the stellar radius, i.e. the lat- 
est when /?2,r = Ri,e- Once mass transfer re- 
sumes the binary's orbital period is P\ < P u . 
In other words: the binary has crossed the pe- 
riod range P\ < P 0I \, < P u as a detached sys- 
tem. And because of lacking accretion lumi- 
nosity, such systems are intrinsically very faint, 
fainter even than pre-CVs, and, therefore, vir- 
tually unobservable. A gap in the period dis- 
tribution can thus arise if a) a sudden drop 
of the AML rate causes CVs to detach, b) if 
that happens to most of the CVs evolving from 
P or b > P u — > Porb < Pu, and if c) the values of 
P u and Pi are practically the same for all sys- 
tems going through a detached phase. 

So far I have not yet addressed the ques- 
tion why the AML rate should drop by a 
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large factor and that also on a sufficiently 
short time scale. The idea behind this propo- 
sition is that effective amplification of mag- 
netic flux via a dynamo and thus efficient 
AML loss via "magnetic braking" is strongly 
tied to the presence of a convective envelope 
and a radiative core in the donor star (e.g. 
ISpruit & Ritter I dl983lf ). Accordingly, it is pro- 
posed that AML via "magnetic braking" de- 
creases rapidly when, as a consequence of on- 
going mass loss, the donor eventually becomes 
fully convective. The questions of whether that 
really happens and whether AML via "mag- 
netic braking" stops completely or only par- 
tially when the donor becomes fully convec- 
tive have remained somewhat controversial to 
this day. Qualitative theoretical arguments in 
favour of the abov e proposition have, how ever, 
been presented by Taa m & Spruit I d 1989b . 

In order for the disrupted magnetic brak- 
ing proposition to work quantitatively the fol- 
lowing requirements must be met: AML above 
the gap must drive mass transfer at a level of 
-M 2 ~ 10 _9 M G yr _I . As a result, the donor be- 
comes fully convective when M conv ~ O.2M 
and P a ~ 3 h . At that moment, as a consequence 
of previous high mass loss, the stellar radius 
is larger by about 30% than in thermal equi- 
librium. With the disappearance of the AML 
from "magnetic braking" the AML loss rate 
drops by a factor of ~ 10 - 20 to essentially the 
value due to gravitational radiation alone. After 
the detached phase which lasts ~ 10 9 yr mass 
transfer resumes with M 2 = M conv ~ O.2M , 
Ri = Ri, e ~ O.2fl and P orb = P x ~ 2 h at a 
level of -M 2 ~ 5 10 _11 M o yr _1 . Explaining the 
gap as a collective phenomenon of CV evolu- 
tion requires furthermore that the majority of 
the donor stars are all of the same type, i.e. 
MS stars, and that AML via "magnetic brak- 
ing" yields similar mass transfer rates in dif- 
ferent systems at the same orbital period. This 
guarantees that P u and Pi are more or less the 
same for all systems and thus the coherence of 
the phenomenon. 

The fact that the period range of the gap 
is not empty already indicates that not all CVs 
follow the above-described evolution strictly. 
There are several reasons for why there may 
be CVs in the gap. The most important ones 



are: 1) a donor mass such that at the end of 
the detached pre-CV evolution the orbital pe - 
riod is 2 h < P mh < 3 h (e.g. iKolbl (119931) . 
iDavis et aTTI (l2008h ~): 2) a donor star which ini- 
tiall y was close t o the terminal age MS (see 
e.g. lRitterl (1 19941) ). or which is the artificially 
evolved remna nt of earlier thermal tim e scale 
mass transfer (ISchenker & King 1120021) : 3) re- 
duced "magnetic braking" because of the pres- 
enc e of a strongly magnetized WD (fo r details 
see ILL Wu & Wickramasinghe I (119941) 1 . 

At the end of CV evolution the donor 
star is a very faint brown dwarf. The WD, in 
turn, with an effective temperature of typically 
< 10 4 K is also very faint. And because the 
mass transfer rate resulting from gravitational 
radiation is very small as well, i.e. -M 2 ^ 
10 _11 M o yr _I , so is the resulting accretion lu- 
minosity. Thus, such CVs are extremely faint 
and inconspicuous objects, and correspond- 
ingly difficult to detect. And though intrinsi- 
cally about 90% o f all CVs are in this late 
phase (iKolb II 19931) so far only one convincing 
candidate beyon d and far from the p eriod mini- 
mum is known f Littlefair et al. Il2006l) . The CV 
graveyard, as this evolutionary branch is some- 
times referred to, is thus largely hidden from 
our view. 
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